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Abstract 



We present a general method to obtain static anisotropic spherically symmetric solutions, satisfying a non- 
local equation of state, from known density profiles. This equation of state describes, at a given point, the 
components of the corresponding energy-momentum tensor not only as a function at that point, but as a 
functional throughout the enclosed configuration. In order to establish the physical aceptability of the pro- 
posed static family of solutions satisfying nonlocal equation of state, we study the consequences imposed by 
the junction and energy conditions for anisotropic fluids in bounded matter distribution. It is shown that a 
general relativistic spherically symmetric bounded distributions of matter, at least for certain regions, could 
satisfy a nonlocal equation of state. 



1 Introduction 



The structure of a relativistic star is believed to be rather complicated (solid crust, superfluid interior, 
different exotic phases transitions). Despite that its bulk properties seem to be computed with reasonable 
accuracy making several simplifying assumptions, the true equation of state that describes the properties of 
matter at densities higher than nuclear (« 10 14 gr/cm. 3 ) is essentially unknown. This is mainly due to our 
inability to verify experimentally the theories that should describe the microphysics of nuclear matter at such 
high densities W, 0, ||, [| . Presently, what is known comes from the experimental insight and extrapolations 
from the ultra high energy accelerators and experimental cosmic physics (see || [fj and references therein) . 
Having this uncertainty in mind, it seems be reasonable to explore what is allowed by the laws of physics, 
in particular, within the framework of the theory of General Relativity and considering spherical and axial 
symmetries. 

Classical continuum theories are based on the assumption that the state of a body is determined entirely 
by the behavior of an arbitrary infinitesimal neighborhood centered at any of its material points. Further- 
more, there is also a premise that any small piece of the material can serve as a representative of the entire 
body in its behavior and, hence the governing balance laws are assumed to be valid for every part of the body, 
no matter how small. Clearly, the influence of the neighborhood on motions of the material points, emerging 
as a result of the interatomic interaction of the rest of the body, is neglected. Moreover, the isolation of an 
arbitrary small part of the body to represent the entire one clearly ignores the effects of the action of the 
applied load at distance. These applied loads are important because their transmissions from one part of 
the body to another, through their common boundaries, affect the motions and hence the state of the body 
at every point. 

The relevance of long-range or nonlocal outcomes on the mechanical properties of materials are well 
known. The main ideas of non local a continuum were introduced during the 1960s and are based on 
considering the stress to be a function of the mean of the strain from a certain representative volume of the 
material centered at that point. From that time there have been many situations of common occurrence 
wherein nonlocal effects seem to dominate the macroscopic behavior of matter. Interesting problems coming 
from a wide variety of areas such as damage and cracking analysis of materials, surface phenomena between 
two liquids or two phases, mechanics of liquid crystals, blood flow, dynamics of colloidal suspensions seem to 
demand this type of nonlocal approach which has made this area very active concerning recent developments 
in material and fluid science and engineering (see J?J and references therein). 

In a recent work it is shown that under particular circumstances a general relativistic spherically 
symmetric anisotropic distribution of matter could satisfy a Nonlocal Equation of State (NLES from now on). 
Some types of dynamic bounded matter configurations having a NLES, with constant gravitational potentials 
at the surface, admit a Conformal Killing Vector and fulfill the energy conditions for anisotropic imperfect 
fluids. We also developed several analytical and numerical models for collapsing radiating anisotropic spheres 
in general relativity. 

The present paper is focussed in a more general family of NLES and in determining the conditions under 
it could represent reasonable bounded matter distribution in General Relativity. 

The static limit of the particular NLES considered in reference (|] can be written as 



P r (r) = p(r) - 1 jT f 2 p(f) df + ^ ; (1) 

where C is an arbitrary integration constant. It is clear in equation (Q) a collective behavior on the physical 
variables p(r) and P r {r) is present. The pressure Pr{r) is not only a function of the energy density, p(r), at 
that point but also its functional throughout the rest of the configuration. Any change in the pressure takes 
into account the effects of the variations of the energy density within an entire volume. 

An additional physical insight of the meaning of the nonlocality for this particular equation of state can 
be gained by considering equation ([!]) re-written as 

*m - «) - 1 Wr + £ , ^ Wr = . , 2) 

Clearly the nonlocal term represents an average over the function p(r) within the volume enclosed by the 
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radius r. Moreover, equation (g) can be easily rearranged as 

P(r) = \p{r) + ~ (p(r) - (p(r)» + ^ = \p(r) + \a p + ^ , (3) 

where we have used the concept of statistical standard deviation a p from the local value of energy density. 
Furthermore, we may write: 



C { V{r) = \p(r) 

P r (r)=V(r) + 2a v{r) + -5 where 1 (4) 

I °v{r) = {\p{r)-k{ P ) r ) = {V{r)-V{r)) 

Therefore, if at a particular point within the distribution the value of the density, p{r), gets very close to 
its average (p{r)) the equation of state of the material becomes similar to the typical radiation dominated 
environment, P r (r) ~ P( r ) = \p{ r )- 

The structure of this contribution is the following: first we give the general conventions and the field 
equations; secondly in Section III the nonlocal equation of state for anisotropic static models is presented; 
following is Section IV where we consider the consequences imposed by the junction and energy conditions; 
and finally the new solutions are shown in the Section V. 



2 The Einstein Field Equations 

To explore the feasibility of nonlocal equations of state for bounded configurations in General Relativity, we 
shall consider a static spherically symmetric anisotropic distribution of matter with an energy-momentum 
represented by T£ = diag (p, —P r , —P±, —Pi)- Here p is the energy density, P r the radial pressure and P± 
the tangential pressure. Although the perfect pascalian fluid assumption (i.e. P r = P±) is supported by 
solid observational and theoretical grounds, an increasing amount of theoretical evidence strongly suggests 
that, for certain density ranges, a variety of very interesting physical phenomena may take place giving rise 
to local anisotropy (see |J and references therein). 

We adopt standard Schwarzschild coordinates (t, r, 9, <fi) where the line element can be written as 



with dtt 2 = d6 2 + sin 2 



d.s 2 = e 2 ^dt 2 
the solid angle. 



s 2A(r) dr a 



r 2 dn 2 



The resulting Einstein equations are: 

8tt P = 

-8TtP r = 
SirP, = 



1 + 


g -2A 




r 


1 


e -2A 




r 


e" 2A 


'A' 

r 



2A'- 
2v' + 



r 



v" + - (z/) 



and 



(5) 

(6) 
(7) 
(8) 



where primes denote differentiation with respect to r. 

Using equations (^) and (||), or equivalently the conservation law T^.^ = 0, we obtain the hydrostatic 
equilibrium equation for anisotropic fluids 



P; = - (p + P r )v' + - (P ± - P r ) 



It can be formally integrated to give 



where the mass function m(r) defined by 



e- 2A = 1 



m(r) 



m(r) — An / pr 2 dr, 



(9) 



(10) 



(11) 
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is the mass inside a sphere of radius r as seen by a distant observer. Finally, fromlQ), (g) and (m) the 
anisotropic Tolman-Oppenheimer-Volkov (TOV) equation pj can be written as 

dP r , s fm + 4irr 3 P r \ 2 , 

—L = -( fi + P r )l— —1 +-P X -Pr ■ (12) 

dr y r(r - 2m) / r 

Obviously, in the isotropic case (P± = P r ) it becomes the usual TOV equation. 



3 A Family of Solutions with a NLES 

In this section we are going to present a family of static solution of the Einstein Equation satisfying a NLES 
Defining the new variables: 



e Mr) =h(T) e 4/5(r) and e 2A(r) = J_ . 

h(r) 

the above metric (^) can be re-written as 



with /i(r) = 1-2 



m (r) 



ds 2 = h (r) e 4 ^dt 2 - -^—dr 2 - r 2 dft 
h{r) 



and the resulting Einstein Equations are: 

l-h-h'r 

07T/9 = 3 

87TP,, 



l-h-h'r Ah (3' 



and 



8ttP x = 



h' + 2hd' 1 



h" + 4 h p" + 6 + 8 h (P' f 



Now, if equation (|l|) is re-stated as 

p - '3P„ +r(p'-P;) = 0, 
we have, from ( p^| ) and by using (|l5|) - (p^|), 

2 
r 

It can be formally integrated yielding 



(h' + 2/i/3') + h" + 20' ti + 2/i/3" = 



1 



C 



c 



Hl > = 2 l »{h) + J— h dr + ^ 

where C and C\ are arbitrary integration constants. 

The corresponding Einstein equations in terms of the metric elements ([l3]) are: 

2m' 



8tt P - 



8ttP = 8irp - 



4 (m - C) 



and 



87rP, 



m"r 2 (r - 2m) + 2r 2 m' (m 1 - 1) - 2m (m - r) - 2C (r + 2m - 3rm') + 4C 2 



r 3 (r — 2m) 

At this point, equation (|20| ) deservers several comments. 

• Firstly, if we set C — C\ = , a particular solution is found, i.e. 



0(r) 



2 I /i 



(13) 

(14) 

(15) 
(16) 
(17) 

(18) 
(19) 

(20) 

(21) 
(22) 
(23) 



(24) 



with C a constant parameter. We have considered this family of solutions in a previous work |p[ and 
some of their geometric properties and several collapsing models of radiating anisotropic spheres were 
presented. 
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• The second comment concerning equation ( p0| ) is the approach we have followed in order to obtain 
static anisotropic solutions having a NLES. It is clear that if the profile of the energy density, p(r), is 
provided, the metric elements h(r) and j3(r) can be calculated through (jll]), (jl^) and (^0|). Therefore we 
can device a consistent method to obtain static solutions having NLES from a known static solutions. 

It is clear that the metric elements describing bounded matter distribution should fulfill the junctions 
conditions and the physical variables coming from the energy momentum tensor are only restricted by some 
elementary criteria of physical acceptability and the hydrostatic equilibrium equation ( |l2] ) . The next section 
is devoted to list those conditions for anisotropic fluids. 



4 Energy and Junctions Conditions 

Most of the spherically symmetric exact solutions found in the literature do not represent physically "real- 
istic" fluids (see for examples, two interesting and complementaries reviews on this subject 11 and jjJJ). 
The condition of what is a realistic fluid is subjective and depends from author to author on which of the 
energy conditions are considered valid. In order to establish the physical aceptability of the proposed static 
family of solutions ( pp| ) satisfying a NLES, we shall study consequences imposed by the junction and 
energy conditions for anisotropic fluids on bounded matter distribution. 
In order to be matched to the exterior Schwarzschild solution, 

ds 2 = (l - 2yj dt 2 - (l - 2yj dr 2 -r 2 d0 2 , (25) 

the interior metric should satisfy the following conditions at the surface of the sphere r = a : 

(3 (a) = p a = 0, m(a) = M and P r (a) = (26) 



thus equation (22) leads to 

C = M — 2na 3 p (a) (27) 
The elementary criteria we use to characterize a physically meaningful anisotropic fluid are: 
1. The pressure must be positive: 

P r > and P± > (28) 



2. The pressure gradient must be negative: 



3. The speed of sound should be subluminal: 



BP, 

°g < (29) 



BP 

v s ^^<l (30) 
op 

4. The trace of the energy-momentum tensor must be positive (Strong Energy Condition): 

p>P r + 2P ± (31) 

5. The density must be larger that the pressure (Dominant Energy Condition): 

p> P r and p > P ± (32) 
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5 The Method and NLES Static Solutions 

In the present section we state a general method to obtain NLES static anisotropic spherically symmetric 
solution from density profiles. 



1. Select a static density profile p(r), from a known static solution. Thus, the mass function can be 
obtained through equation (|ll| ) and by the junction condition (|26| ) fulfills the continuity of h(a) giving 
the expression for the total mass m(a) = M . 

2. Following, the metric coefficient /3(r) can be found by using equation (|2 



3. Finally, Einstein equations (J22J) through (|23j) provide the expressions for the radial and tangential 
pressures, P r and P+, respectively. The integration constants are obtained as consequences of the 
junction conditions at the boundary surface, r — a, i.e. /3(a) = 0, 771(a) = M and P r (a) = 0. 

4. Once we have a candidate for a bounded matter distribution having a NLES we should explore where, 
within the configuration, all the above conditions (equations (28) through (p3)) are valid. 



In order to illustrate the above procedure we shall work out seven static solutions: Tolman VI, and 
Tolman- V Jig)] , an static solution proposed by B. W. Stewart [14|], three solutions proposed by M. C. 
Durgapal ]15| and the M. K. Gokhroo and A. L. Mehra solution (XC J , (see Appendix). Tolman and Durgapal 
solutions are borrowed from isotropic static distribution of matter but Stewart and Gokhroo & Mehra 
solutions describe static bounded anisotropic ones. Stewart solution outlines a conformally flat anisotropic 
distribution which is more stable than its isotropic counterpart. 

The Gokhroo and Mehra solution corresponds to an anisotropic fluid with variable density and under 
some circumstances |l^] , represents densities and pressures given rise to an equation of state similar to the 
Bethe-Borner-Sato newtonian equation for nuclear matter [[l], ^| . 

The expressions for the starting energy density can be summarized as follows: 

Singular models: 



No singular models: 



Equation of State 



Tolman V 



Tolman VI 



Stewart 



Density p 



1 

8 71 



3 
7r 2 



+ 



3R 2 



' jcr 

~ (e 2Kr -l)(e ikr - +8Kre 2Kr_ 1 ) 



8nr 2 



Equation of State 



Density p 



Durgapal (n=l) 



l-3K-3Kx 
(1+2 x) 



2(1 + Kx) 
(1+2 x) 2 



Durgapal (n=2) 



C K(3+5~xJ ~ 

87r (1+3 x)i 



Durgapal (n=3) 



c 



87r(l+a:) 2 



|(3 



x _ 3K(l+3x) 
(1+4 a:) 3 



Gokhroo &: Mehra 



The parameters that characterize bounded configurations: mass M, M/a, boundary redshift z a , surface 
density p a and central density p c the are the shown in the following table. Let us consider stars with radius 
of 10 Km. 



Equation of State 


M/a 


M (Mq) 


Za 


p a x 10 14 (gr.cm a ) 


TolmanV 


0.25 


1.69 


0.4 


3.57 


TolmanVI 


0.25 


1.70 


0.4 


2.68 


Stewart 


0.32 


2.15 


0.6 


6.80 



and for the models with central density: 
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Equation of State 


M/a 


M (Mn) 
V W } 




p a x 10 14 (gr.cm 3 ) 


p c xl0 ib (gr.cm a ) 


Durgapal (n=l) 


0.25 


1.69 


0.4 


5.36 


2.00 


Durgapal (n=2) 


0.33 


2.26 


0.7 


7.15 


2.70 


Durgapal (n=3) 


0.38 


2.54 


1.0 


8.04 


2.41 


Gokhroo & Mehra 


0.20 


1.35 


0.3 


4.29 


0.96 



Within the above configurations the range of validity for the NLES are 



Equation of State 


r (Km) 


Tolman V 


7.8 < r < 10.0 


Tolman VI 


7.5 < r < 10.0 


Stewart 


0.7 < r < 10.0 


Durgapal (n= 1 ) 


0.0 < r < 10.0 


Durgapal(n=2) 


0.0 < r < 10.0 


Durgapal(n=3) 


0.0 < r < 10.0 


Gokhroo & Mehra 


1.1 < r < 8.2 
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7 Appendix 



7.1 Tolman V-like solution 

We shall now take the solution V of Tolman as providing the first of the examples of a sphere of fluid 
surrounded by empty space. The density is given by: 



1 

s7 



" 3 
Tr 2 " 



10 
3~R2 



(*)*' 



(33) 



^From equations (jll]) and ([!]) we obtain the expressions for the mass and radial pressure, 



+ 

7r 2 R 2 \R 



and P r 



An 



respectively. 

^From the boundary conditions, C± and R are found to be 



84 



•-»(s) ! 



and R = 



3 J_/_^ 

1^ + R?\R 



_ 3 10 

7?a~ 



27rr 3 



(34) 



(14M - 3a) 



(35) 



The following table contains the ranks of validity for the different criteria indicated in the section IV 





p>0 


Pr > 


P± > o 


^<0 

Or 


m. < i 


p > P r + 2 P ± 


P > Pr 


P>P± 


r e 


(0, 10.0] 


(0, 10.0] 


(0, 0.3] 

u 

[7.8, 10.0] 


(0, 10.0] 


[5.8, 10.0] 


[3.4, 10.0] 


[3.6, 10.0] 


[2.0, 10.0] 



7.2 Tolman Vl-like solution 

In this case, the model to be studied is the solution VI of Tolman. We recall that the equation of state of 
this model, for large p, approaches that for a highly compressed Fermi gas. The density matter is 



56 7r r 2 

Respectively, the mass and pressure are found to be 



P=J^ (36) 



m = T-crr and P r = — p + - . (37) 
14 2-Kr 6 



The constant C\ and the boundary conditions lead to 



3 , 14M 

C = —aa and a = — — 38 
28 3 a 



The energy conditions yield to 



r < 10.0 


p>0 


Pr > 


P± > 


^f<0 

or 


dp - 


p> Pr+2P ± 


P> Pr 


P>P± 


r e 


(0, 10.0] 


(0, 10.0] 


(0, 5.0] 


(0, 10.0] 


[7.5, 10.0] 


[5.0, 10.0] 


[5.0, 10.0] 


[3.0, 10.0] 
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7.3 Stewart Solution 

B.W. Stewart examines several anisotropic, conformally flat, internal solutions. He presents four examples 
of mass distributions, we select the example 2 from that paper E3] 



P : 



i 



,2Kr 



1) 



8Kre 2Kr - l) 



nr z (e 2Kr + 1) 

The expressions for the mass distribution and pressure are 



(39) 



r l e 
m = — 



2 AY 



- 1 



and 



Pr =P 



1 



JlKr 



1 



Airr 2 \ e 2Kr + 1 



Again, the constant C\ is obtained to make P r (a) — 0, i.e. 

_ a e 2Ka (e 3Ka - 8Kae 2Ka - e 2Ka + 8Ka - l) 
C ~ 4 e 2Ka (e 3Ka + 3e 2Ka + 3) + 1 

while the constant K is found from the boundary conditions as 



2 7rr 3 



(40) 



(41) 



K = — In 

2a 



1+(2 M )5 



'2M\ 2 



(42) 



The energy conditions lead 



r < 10.0 


p> 


Pr > 


Pi. > 


^<0 

or 




p > P r + 2 Pjl 


P > P- 


P> Pi 


r e 


(0, 10.0] 


(0, 10.0] 


(0.5, 10) 


(0, 10.0] 


[1.9, 10.0] 


[0.6, 10.0] 


(0.5, 10.0) 


(0, 10.0] 



7.4 Gokhroo and Mehra solution 

Now, the Gokhroo and Mehra |l(J solution is presented. The mass density is 

„2 " 



p = a 

the mass and the radial pressure, can be written as 

.3 



1 - K- 



a r 
~6 



1 



3Kr 2 



3 1 
and P r = -/?- — 

5 Z7T 



J7_ _C 

15 ~ r 3 



and the boundary conditions lead to 



C = ^![9X-5] and o= J° M - 
60 L J a 3 (5-3X) 



The physical variables are valid in the following range 



(43) 



(44) 



(45) 



r < 10.0 


P >0 


Pr > 


Pl > 


^<0 

or 


Op - 


p > P r + 2 P ± 


P > Pr 


p>p± 


r e 


(0, 10.0] 


(0, 10.0] 


(0, 8.2] 


(0, 10.0] 


[1.1,10.0] 


(0, 10.0] 


(0, 10.0] 


(0, 10.0] 



7.5 Durgapal Solution 

The paper of M. C. Durgapal jl5| report several solutions. The autor assume that the value of the metric 
coefficient is given by: e u — A (l + Cr 2 ) , where A and C are constant and n a parameter. Durgapal shows 
that Einstein's equations can be solved explicitly for different values from n. 
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• For n = 1 

The first solution are identical to Tolman IV solution. The expression for the density is given by 

C 



9 : 



8tt 



1 - 3K - ZKx 2(1 + Kx) 
(l + 2a;) + (l + 2x) 2 



(46) 



Again, the expressions for mass and radial pressure are obtained from the equation (|l|): 

xi K(l + x)-l , C K(l+x)-l C 

m = ; — — and P r = p 



where 



2C5 (l + 2x) 



x = C r: xi — C a = 



M 



o-3M' 



and K = 



An (l + 2x) 27r(§) ! 
C[axi -2M(l + 2x x ) 



x\ (1 + xi) 



(47) 



(48) 



The constant K is obtained for boundary conditions. The pressure in the surface (P r {a) = 0) leads to 
that: 



C 



axi K(l + x 1 + 2x{) +2a;i-l 
4 l + 4ari(l + a;i) 



(49) 



The validity range are: 



r < 10.0 


P>0 


P r > 


P± > o 


^<0 

or 


dp - 1 


p>P r + 2P x 


p > p r 


p>p± 


r € 


(0, 10.0] 


(0, 10.0] 


(0, 10.0] 


(0, 10.0] 


(0, 10.0] 


(0, 10.0] 


(0, 10.0] 


(0, 10.0] 



For n = 2 

This solution is identical to that obtained by Kuchowicz, Adler and Adams and Cohen, according to 
the author. The density is given by 



P 



C K(3 + 5x) 



8tt 



(l + 3aO J 



The mass and radial pressure, are : 



1 Kx? 



2C2 (i + 3x)i 



and P r = p 



1 CK 



4tt 



(l + 3a;) 3 2 7 r(§) i 



where 



x = C r; x\ — C a — 



M 



2M a 

and K = (l + 3xi) 3 

ax\ 



2a- 5M' 

The constant K is obtained for boundary conditions and C by (P r (a) = 0) 

_ CKa 3 1 - X! 
4 (l + 3a;i)* 

The regions where the physical variables are valid 



(50) 



(51) 



(52) 



(53) 



r < 10.0 


p>0 


Pr > 


P± > 


^<0 

or 


^<1 


p > P r + 2 Pi 




p>p± 


r € 


(0, 10.0] 


(0, 10.0] 


(0, 10.0] 


(0, 10.0] 


(0, 10.0] 


(0, 10.0] 


(0, 10.0] 


(0, 10.0] 



For n = 3 

In the third solution studied by Durgapal the density has the following form: 

C 



8tt(1 + xY 



| (3+I) _^(i±l£) 

2 (l + 4x) 5 



(54) 
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This expression for the density is exactly the same expression for the densited reported in Durgapal 
and Banncrji for the case with K = 

Again, the mass and the radial pressure are 



m = — - 



2K-3(l + 4x)* 
4Cs (1 + x) (1 + 4x)5 



and 



P + 



C 



where 



x = C r; X\ — C a 2 



M 



and K 



C 2K- 3(1 + Ax) 1 
(1 + 4^)5 



3a-7M' 2a X! 

Again, tc constant K is obtained for boundary conditions and C by: (P r (a) = 0) 

x\ 2K(l- Xl -&x 2 ) -3(1 + 4a;i)5 (1 + 3x1 -4a;?) 



8Ci 



(l + 4xi) 2 (l + 2a;i +xf) 



The validity ranks are 



(55) 



[4M(l+xi)-3oxi] (56) 



(57) 



r < 10.0 


p>0 


Pr > 


P± > o 


^<0 


t)K < j 
a P - 1 


P > P r + 2P ± 


P> Pr 


P>P± 


r e 


(0, 10.0] 


(0, 10.0] 


(0, 10.0] 


(0, 10.0] 


(0, 10.0] 


(0, 10.0] 


(0, 10.0] 


(0, 10.0] 
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